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Abstract. We calculate the normal metal-s+g-wave superconductor tunnelling 
spectrum for various junction orientations and for two forms of the superconducting 
gap, one which allows for point nodes and the other which allows for line nodes. 
For a junction oriented with its normal parallel to the ab plane of the tetragonal 
superconductor, we find that the tunnelling spectrum is strongly dependent on 
orientation in the plane. The spectrum contains two peaks at energies equivalent 
to the magnitudes of the gap function in the direction parallel to the interface normal 
and in the direction making a 7r/4 angle with the normal. These two peaks appear in 
both superconductors with point nodes and line nodes, but are more prominent in the 
latter. For the tunnelling along the c axis, we find a sharp peak at the gap maximum in 
the conductance spectrum of the superconductor with line nodes, whereas with point 
nodes we find a peak occurring at the value of the gap function along the c axis. We 
discuss the relevance of our result to borocarbide systems. 
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1. Introduction 

Nonmagnetic rare-earth borocarbides, such as YN^E^C and LuN^E^C, are among the 
materials which exhibit unconventional superconductivity. There is strong evidence 
in these materials indicating that their superconducting gap function is highly 
anisotropic [TJ E] and that there exist low-lying excitations in the superconducting 
state [3 |U 13 El • The presence of these low energy excitations implies that the gap 
function has nodes (or deep minima) on the Fermi surface. The location of these nodes 
and their nature, i.e., whether they are point nodes or line nodes is still unclear. The 
thermal conductivity in a magnetic field at low temperatures [JJ and the dependence 
of the specific heat on the magnetic field suggest the existence of line nodes 
like those in the cuprates and UPt3. However, the recent measurements of the c-axis 
thermal conductivity in a rotational magnetic field along the ah plane were interpreted as 
evidence for a gap function with point nodes along [100] and [010] directions [2|. Further 
studies which employ different experimental techniques may be necessary to determine 
the detailed structure of the superconducting gap function in momentum space. Here 
we suggest that directional tunnelling spectroscopy may be useful for this purpose. 

In general, the tunnelling conductance spectrum of an anisotropic superconductor 
is strongly dependent on the crystal orientation with respect to the interface plane. 
In the case of a c?-wave superconductor with vertical line nodes, it has been shown 
that features in the conductance spectrum occur at voltages which depend strongly 
on crystal orientation. These voltages correspond to values of the gap function at 
particular points on the Fermi surface For instance, in the case of an ah plane 
tunnelling junction, if the surface orientation of the superconductor is not [100] or [010], 
the conductance spectrum contains a peak at the voltage corresponding to the value 
of the gap function in the direction parallel to the surface normal. The observation of 
this feature would, in principle, allow the directional normal metal-superconductor (NS) 
tunnelling spectroscopy to map out the magnitude of the superconducting gap function. 

In this paper, we calculate the NS tunnelling spectra of anisotropic s-wave 
superconductors with two forms of gap functions, which have been suggested as possible 
candidates for the gap in borocarbides jTOJ ITT] . The candidates are allowed by the 
symmetry of the borocarbide crystal structure and consistent with experimental results. 
It is apparent that the gap function has four-fold symmetry and does not change sign but 
there is evidence for deep minima along [100] and [010] directions [JJEIIHJIUEJIIIJIIIIB]- 
The simplest suitable candidates are of the s + g-wave form. In our calculation, we find, 
similar to a <i-wave case, that features appear in the spectra at voltages corresponding 
to values of the gap at particular points on the Fermi surface depending on surface 
orientation. Some of these features are characteristic of the general s + g-wave gap 
function and are thus common to both candidates, while one particular feature occurs 
only for the candidate with point nodes. Thus, if these features can be observed 
experimentally, it may be possible to determine which form of the gap, if either, correctly 
describes that of the borocarbides. 
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The two candidate gap functions are both s + g-wave and are given by equation (|lqj) 
and (flTijl respectively. 



where 9 and <fr are the polar angles in spherical coordinates, and A s and A g are the s and 
g components of the gap function respectively. For an elliptical Fermi surface, the case of 
line (point) nodes occurs when A s = A g in equation ([Taj) (equation (flftjl ). We study both 
ah plane and c-axis tunnelling spectra for both forms of the gap function by using the 
so-called Blonder-Tinkham-Klapwijk (BTK) scattering formalism ^2] and a continuous 
model to describe the electronic structure of the normal metal and superconductor. 
Although detailed features of the conductance depend on the shape of the Fermi surface, 
it is sufficient to use the continuous model which gives an elliptical Fermi surface to study 
the positions of the main features in the conductance spectra 0. 

In the case of tunnelling into the ah plane of a s + g-wave superconductor, we show 
in this paper that there are three main features occurring at voltages associated with 
the difference of the s and g component of the gap function and the values of the gap 
function in two particular directions with respect to the interface normal, one in the 
direction parallel to the surface normal and the other in the direction making a 7r/4 
angle with the surface normal. These two features are more prominent for the s + g- 
wave superconductor with line nodes than for the superconductor with point nodes. In 
the case of c-axis tunnelling junctions, there is one prominent feature at the maximum 
gap in the conductance spectrum of the superconductor with line nodes, whereas for the 
superconductor with point nodes the conductance spectrum contains two features, one 
at the maximum gap and the other at the value of the gap function along the c axis. 

In Section El we describe the calculation method and the assumptions used 
throughout this paper. Then, we provide the detailed results and discussion of all 
the cases of interest at zero temperature in Section El In addition to the cases in which 
the nodes exist, we also consider cases in which the superconducting gap has small but 
nonzero gap minima. We show all the results in both the Andreev limit (low barrier) 
and the tunnelling limit (high barrier). 

2. Assumptions and method of calculation 

As in reference ^2], we represent the NS junction with an infinite system, the left half 
of which is a normal metal and the right half of which is a superconductor (see figure 
The insulating barrier is represented by a delta function potential with strength H. For 
the ab plane tunnelling junctions, the interface normal vector lies somewhere in the ah 
plane, and for the c-axis tunnelling junction the interface normal is parallel to the c 
axis. 

We take the normal metal to be cubic, and take the superconductor to be tetragonal 
to describe the crystal structure of the borocarbides. In our calculation, we ignore both 



A fcjl = A s - A 9 cos4</> 
Afc ; 2 = A s — A g sin 2 9 cos 40 
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Figure 1. The normal metal-superconductor (NS) junction is represented by an 
infinite system as shown in this picture. The normal metal fills the x < region 
and the superconductor fills the x > region. The insulator layer is represented 
by a delta function of height H in units of energy per length. The gap function is 
taken to be zero in the normal metal and to be nonzero and independent of x in the 
superconductor. 



the suppression of the gap function near the NS interface and the proximity effect for 
simplicity. The gap function A k is taken to be as in either equation (JTaJl or (JT5J). 

The Bogoliubov-de Gennes equations that describe the excitations of the system 

are 



U(r) = EU(r) 



6 p + H5(x)-n A k B(x) 

A k Q(x) -6 P - H6(x) + n 

where \i is the chemical potential, 0(x) is the Heaviside step function, 
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m (effective mass in the normal metal), x < 

m a b(ab plane effective mass in the superconductor), x > 

m, x < 

m c (c-axis effective mass in the superconductor), x > 



and U(r) is a two-component function: 
U(r) 



u(r) 




Uk 


v(f) 







ik-r 



(3) 



After substituting U(r) from equation (JH} into equation^, we obtain the bulk excitation 
energies for the normal metal and superconductor respectively as 

i<7 



E{k) 



e k + a 2 



(4a) 
(46) 
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Figure 2. The plots of two excitation energies as a function of q x or k x , the component 
parallel to the interface normal, at a particular value of fcii, which is (k y ,k z ). The 
plot on the left is for the excitation energy of the normal metal and the plot on the 
right is for the excitation energy of the superconductor. At the same energy, there 
can be 4 propagating excitations for each side. However, for an incoming electron 
from the normal side, the wave function of the normal metal is a linear combination 
of only 3 excitations represented by the open circles, and the wave function of the 
superconductor is the sum of 2 outgoing excitations also represented by the open 
circles. 



where for the normal metal (equation (j4a)) )the plus and minus signs are for electron and 
hole excitations respectively, and for the superconductor (equation (|46l)) 

^2 / i.2 _l 12 U 2 

~2 



K + K , K 



-/!. 



(5) 



m ah m c 

Figure 121 shows a plot of the excitation energy of the normal metal (superconductor) 
as a function of q x (k x ), the component along the interface normal, at a particular 
q[\ — k\\ — (k y , k z ), the component perpendicular to the interface normal. 
The amplitudes of the excitations, Uk and Vk, of the normal metal are 
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for electrons 



for holes 



(6) 



whereas those of the superconductor are 

Uk 1 



y/\E + £ k \* + A h 



E + £ k 



(7) 



The wave function of each side is a linear combination of all the appropriate excitations 
of the same energy and the momentum that has the same component perpendicular 
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to the interface normal. For ah plane tunnelling, the wave functions of both sides are 
therefore 



U N (r) 



where q* 1 and fc ± (see figure EJ) satisfy 
hq ± = 



1 
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(96) 
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and a, 6, c, and d are the Andreev reflection, the normal reflection, the same-branched 
transmission, and the cross-branched transmission amplitudes respectively. For the c- 
axis tunnelling, we can obtain the wave function of each side in a similar way. 

Normally, the range of the energy E relevant to the NS tunnelling experiments is 
of order meV whereas the Fermi energy is of order eV. Therefore, we use the following 
approximation for q ± and k^: 
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(10a) 
(106) 



where qF is the magnitude of the Fermi wave vector of the normal metal, and /c^, m* are 
the parameters that satisfy h 2 kp/(2m*) = Ep,s, the Fermi energy of the superconductor. 
Using the conservation of the momentum parallel to the surface, we have the following 
relationship between the polar angles in spherical coordinates: 



qF sin 9n sin (pN = kp 
qF cos 6n = kp 



m* 



sin Or sin < 




COS0Q- 



'11a) 



[lib) 



We obtain all the amplitudes a, 6, c, and d by applying the following matching 
conditions at the interface: 
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where Z = mH/(h 2 kF). Note that both matching conditions are for the a6 plane 
tunnelling junction. For the c-axis tunnelling junction, the first condition remains the 
same, but we have to replace with m c in the second condition. 



Normal metals +g -wave superconductor tunnelling spectroscopy 7 

In the BTK formalism, the Andreev reflection, the normal reflection, and the two 
transmission probabilities are used to obtain the current across the junction. All the 
reflection and transmission probabilities are obtained from 

A = \a\ 2 (% J (13a) 



|2 



B = \b\ 2 (136) 
C = \c\ 2 (\u k+ \ 2 - \v k+ \ 2 ) (^j (13c) 

D = \d\ 2 (\u-„-\ 2 - \v. k -\ 2 ) (£) (13d) 



and satisfy: A + B + C + D — 1, i.e., the number of particles is conserved. 

On the normal metal side, we find the current across the junction as a function of 
an applied voltage is 

eQ f 

lNs{V) = J2^J d * VqAl + m ~ Bmf ^ ~ eV) ~ /W] (14) 

where Q is volume, v qx is the x component of the group velocity of the incoming electron, 
and f(E) is the Fermi-Dirac distribution function. The conductance of the ab plane 
junction at zero temperature is 

J Tab 

nab /t/\ _ al NS 



dV 
me^VlV 



>N 



d9 N sin 2 6 N cos (j) N [l + A(V,(J) N ,8 N )-B(V,(j) N ,6 N )] (15) 



A7T 2 h 2 

Similarly, the conductance of the c-axis tunnelling junction is 
me 3 QV f f 

G C NS (V) = J #jv J de N ^6 N cmO N [l + A(y,<j> N ,6 N )-B(y,<j> N ,e N )] 

The limits of both integrals can be found by considering equation ((TT^) and (jllfejl . 



3. Results and discussion 

We plot the normalized conductance as a function of applied voltage of both ab plane 
and c-axis tunnelling. We define the normalized conductance as the conductance of 
the junction normalized by its value at a high voltage, i.e., eV 2> A max , the maximum 
magnitude of the gap function. Using both forms of the gap function in equation (jla|) 
and (|T6j) . we consider 3 cases: (1) A s = A g , (2) A s = 0.9A g , and (3) A s = 1.1A 9 . 
These choices of the parameter A s /A g span the range allowed by the results of thermal 
conductivity measurements in borocarbides which indicate that the ratio of the gap 
maximum to the gap minimum is at least ten pQ. All the results are obtain for zero 
temperature. 
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Figure 3. The geometry of the ab plane tunnelling NS junction. The angle between 
the x axis and the a axis, a, defines the orientation of the junction. 

3.1. ab plane tunnelling 

Figure El shows the diagram of the junction that has the interface normal on the ab 
plane of the superconductor. We specify the orientation in the ab plane with a, the 
angle between the interface normal and the a axis of the superconductor. The gap 
function is a function of a, 

A fc ±,i = A s - A 9 cos4(0 5 ^a) (17a) 

A fc±i2 = A s - A 9 sin 2 # s cos4(0sTa) (176) 

In the Andreev limit (small Z), the conductance spectrum depends very little on 
the junction orientation. As shown in figure |U for A^i and in figure El for A^, the 
conductance spectra in all cases have the inverted gap structure. Note that in the case 
of Afc >2 in figure El there is a slight kink at the voltage associated with the value of the 
gap function along the c axis (eV = A s ). Also note that when A s ^ A g , there is a 
feature occurring at eV = |A S — A s | as marked by the arrows in figure 0] (b),(c) and 
figure El (b),(c). 

In the tunnelling limit (large Z), the shape of the conductance spectrum changes 
with the interface orientation. First, consider the case where A s = A g , there are two 
distinct peaks which are more pronounced for the superconductor with the gap function 
Afc ! than for the superconductor with the gap function Afc 2 (compare figure El (a) and 
figure (a)). The two peaks occur at voltages corresponding to A^^s = n/2,(ps = 
0,a) (marked by the filled arrows) and Ak^(6s = vr/2,05 = 7r/4,a) (marked by the 
hollowed arrows), where i is either 1 or 2. Notice that when a = tt/8, there is only 
one peak due to the fact that Ak^(6s = tt/2, 4>s = 0, a = vr/8) = Ak^(6s = ir/2,<f)s = 
7r/4, a = 7r/8). 

In the case of A s = 0.9A g , A^j can be both positive and negative. Consequently, 
for the orientations with a ^ 0,7r/4, in addition to the two peaks at the positions like 
those in the case of A s = A g , there exists a peak at zero voltage (see figure El (b) and [3 
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Figure 4. ab plane tunnelling junctions. The plots of the normalized conductance 
spectra when Z = for 3 surface orientations in 3 cases of the gap function 
A fe±)1 = A s - A 3 cos4(>s t a): (a) A s = A g , (b) A s = 0.9A S , and (c) A s = 1.1 A g . 
The arrows in (b) and (c) indicate the feature occurring at voltage corresponding to 
| A s — A s |. The insets in (b) and (c) are the enlargements of the conductance plots 
near zero voltage. 
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Figure 5. a6 plane tunnelling junctions. The plots of the normalized conductance 
spectra when Z = for 3 surface orientations in 3 cases of the gap function 
A fc ±, 2 = A s - A g sin 2 S cos 4(<j> s t a): (a) A s = A g , (b) A, = 0.9A 3 , and (c) 
A s = l.lA g . The arrows in (b) and (c) indicate the feature occurring at voltage 
corresponding to \A S — A g \. The insets in (b) and (c) are the enlargements of the 
conductance plots near zero voltage. 
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Figure 6. ab plane tunnelling junctions. The plots of the normalized conductance 
spectra when Z = 3 for 3 interface orientations in 3 cases of the gap function 
Ak± fX = A s - A g cosi{(f>s T a): (a) A, = A S) (b) A s = 0.9A 3 , and (c) A, = 1.1A 9 . 
The hollowed arrows indicate the features occurring at voltage corresponding to 
Afc,i(0s = n/2,<j>s = 7r/4, a) and the filled arrows mark the features at voltage 
corresponding to Ak,i(0s = 7r /2, 4>s — 0, a). 
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Figure 7. a6 plane tunnelling junctions. The plots of the normalized conductance 
spectra when Z = for 3 interface orientations in 3 cases of the gap function 
A fc±i2 = A s - A s sin 2 6»sCOs4(< ? !) S t a): (a) A s = A 9 , (b) A, = 0.9A g , and (c) 
A s = l.lA g . The hollowed arrows indicate the features at voltage corresponding 
to A/- ]2 (#s = 7r/2,05 = 7r/4, a) and the filled arrows mark the features at voltage 
corresponding to A^^s = tt/2, cj)s = 0, a). 




Figure 8. This picture shows a sketch of the gap function in the momentum space, 
when a 7^ 0,7r/4. The dashed circle is the contour of the Fermi surface projected on 
the k x k y plane. The solid curve represents the gap function. The two dashed lines 
are the <ps — 7r/4 and (f>$ = 3ir/4 lines. The solid horizontal line above k y = is the 
line of constant k y at <ps = ir/4 on the Fermi surface. This line cute the Fermi surface 
at 2 points on the plane. Note that d\ = di (the magnitudes of the gap of the two 
excitations when <j)g = 0) and d% = d^ (the magnitudes of the gap the two excitations 
when 4>s = tt/4). 



(b)). The existence of this zero-bias conductance peak, which also occurs in the ci-wave 
case, is a signature of the presence of a sign change of the gap function jTHJ [T3j. 

When A s = 1.1A 9 , there is a finite gap minimum. The conductance is very small 
for voltages less than the gap minimum, and the two peaks at A)~^(9s = 7r/2, — 0, a) 
and Ak t i(9s = tt/2, (fis = tt/4, a) are still present. 

In summary, the conductance spectrum in every case where a ^ 0,7r/4 contains 
two peaks at Ak,i(9s = 7r/2,(fis = 0,a) and Ak t i(9s = vr/2,05 = 7r/4, a). These 
two angles, (fis = 0,7r/4, are special because the magnitudes of the gap function of 
the two transmitted superconducting excitations with the k y corresponding to these 
angles are the same, i.e., A k + = A_ k -, as shown pictorially in figure |S1 The feature at 
^■k,i(9s = 7r/2, 05 = 0, a) has also been shown to occur in a ci-wave superconductor in 
a continuous model [EL However, the feature at = 7r/2,05 = 7r/4, a) is unique 

to an anisotropic s-wave superconductor in this model. In a d-wave superconductor the 
values of A k +(9s = n/2, cj)s = 7r/4, a) and A_ fe -(6>5 = n/2,(j)s = 7r/4, a) always have 
opposite signs and thus excitations having these momenta contribute to the zero energy 
surface bound state [T^] instead of giving rise to such a feature. We note that the two 
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peaks at Ak t i(0s = tt/2, <ps — 0, a) and A^^Os = 7r/2, <ps = tt/4, a) are more prominent 
for the superconductor with the gap function A*.^ than for the superconductor with 
the gap function Afc j2 . Given only this quantitative distinction between the ah plane 
tunnelling results for the two forms of the gap considered, it would likely be difficult 
to determine from real tunnelling data which form is present. The ah plane tunnelling 
spectrum would enable a determination of the magnitude of the gap function in an 
arbitrary direction in the plane for either of the two forms of the gap. 

3.2. c-axis tunnelling 

Because in the continuous model the Fermi surface of a tetragonal crystal is invariant 
under rotation around the c axis of the crystal, the c-axis tunnelling spectroscopy is 
independent of the rotation around the c axis. In the Andreev limit, when |A S — A ff | < 
eV < A max , the conductance curve of the superconductor with the gap function A fc l is 
upward and decreasing smoothly (see figure El (a)). On the contrary, the conductance 
curve of the superconductor with the gap function Afc 2 is downward when |A S — A s | < 
eV < A k>2 (6 s = 0) (see figure M (a)). 

In the tunnelling limit, the spectrum of the superconductor with the gap function 
Ak t i contains a peak at the gap maximum and a feature at | A s — A g \ (see figure El (b)). 
The spectrum of the superconductor with the gap function A^ contains a feature at 
the gap maximum, and a sharp peak at A/,^{0s — 0) (see figure ITU1 (b)). The occurrence 
of the sharp peak at different positions in the conductance spectrum makes it possible 
to distinguish the tunnelling spectrum of the line-node gap from the point-node gap 
form using c-axis tunnelling data. 

4. Conclusions 

We have studied the c-axis and ab plane tunnelling spectroscopy of s + g-wave 
superconductors. The observation of the predicted features in tunnelling measurements 
made for various junction orientations would, in principle, provide a way to study the 
detailed momentum dependence of an s+g-wave superconducting gap. In borocarbides, 
the presence of an s + g-wave gap function having either line or point nodes has been 
suggested. Directional tunnelling spectroscopy would help to determine whether either 
form of the superconducting gap is correct. 

We can distinguish the c-axis tunnelling spectra of a superconductor with line nodes 
from the spectra of a superconductor with point nodes. In the tunnelling limit, the c- 
axis tunnelling conductance spectrum of a line-node superconductor contains a sharp 
peak at the gap maximum, whereas in the spectrum of a point-node superconductor a 
sharp peak occurs at the value of the gap function along the c axis. 

The ab plane tunnelling spectra may be used to map out the magnitude of the 
gap function in the plane. The conductance spectrum of ab plane tunnelling junction is 
strongly dependent on the junction orientation. Two features in the spectrum appear 
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(a)Z 



A=A g 
A=0.9A g 
A=1.1A. 





Figure 9. c-axis tunnelling junctions. The plots of the normalized conductance 
spectra for (a) Z — and (b) Z = 3 in case of the gap A^i = A s — A g cos4(<^s T <*). 
The arrows indicate a feature occurring at |A S — A g \. 



at energies equal to the magnitude of the superconducting gap in two momentum 
directions: one is the direction parallel to the interface normal and the other is the 
direction making a ir/4 angle with the interface normal. It is worth noting that these 
features are more prominent in the spectrum of the superconductor with line nodes than 
in the spectrum of the superconductor with point nodes, although this subtle distinction 
may not be observable especially at finite temperatures. 
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(a) Z = 




(b) Z = 3 



eV/A 
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Figure 10. c-axis tunnelling junctions. The plots of the normalized conductance 
spectra for (a) Z = and (b) Z = 3 in case of the gap At 2 = A s — A s sin 2 9 s cos4(</>g=F 
a). 
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